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We present the results of a perturbative calculation to match the axial and vector currents for 
semileptonic B D'*' decays in lattice NRQCD to the continuum MS scheme. The matching is 
performed to 0{asjY) in Feynman gauge and in the on-shell renormalisation scheme. 

The spatial and temporal components of the currents renormalise differently; to this order the 
matching involves a straightforward renormalisation for the Vb and Ak currents, and a rank two and 
four mixing matrix for the Ao and Vk currents respectively. The resultant one loop corrections are 
of 0(5%), boding well for the accuracy of forthcoming simulations. 



I. INTRODUCTION 

Forthcoming B physics experiments such as BaBar, Belle, CLEG III and LHC-b will greatly 
improve the accuracy with which we know the rates of different decay modes of the B meson, allowing 
improved determinations of the CKM matrix elements involving the 6 quark. The exclusive decay 
B — + DW gives access to Vet- For example, at BaBar the decay amplitude at zero recoil will 
be determined with a statistical error of less than 2% and with similar level of systematic errors, 
a factor of three improvement over CLEG and LEP data. In order to extract Vcb we must divide 
out the form factor at zero recoil, introducing associated theoretical uncertainties. The form factors 
for B ^ D decays are calculable within lattice QCD and it will be highly challenging for the 
community to produce model independent predictions of the form factors with sufficient statistical 
and systematic precision that theoretical uncertainty matches the improved experimental data. 

The decay amplitudes can be factorised into hadronic form factors and perturbative electroweak 
terms as follows: 
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These are often reparametrised in terms of form factors in the heavy quark effective theory, 
{D{p')\V''\B{p)) = VMbMd {h+{v ■ v'){v + v'y + h-{v- v'){v - v'Y] 
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The weak (vector and axial) currents in a lattice-regulated theory differ from the continuum MS 
currents via their different ultra-violet cut-offs. A correction between the results in the two regu- 
larisation schemes is therefore required in order to use simulated lattice hadronic currents in the 



are a number of recent reviews of the status of lattice calculations ^ 
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more commonly quoted continuum MS scheme. Lattice QCD and MS have the same infra-red be- 
haviour, so that these correction factors should be dominated by the differing ultra-violet behaviour 
and one expects that these corrections are calculable within perturbation theory. Knowledge of 
these renormalisation constants enables the various elements of the Cabibbo-Kobayashi-Maskawa 
weak mixing matrix to be probed using experimental decay rates and lattice calculations of hadronic 
matrix elements in conjunction with standard electroweak perturbation theory. 

This paper aims to make possible the accurate determination of B — > D matrix elements by 
showing how to construct the MS currents when using the 0{jj) lattice NRQCD action to simulate 
both the b and c quarks. We perturbatively match the corrections to the weak currents between 
on-shell quark states in lattice NRQCD to the same current in continuum MS with fully relativistic 
quarks, working consistently to 0{jj) throughout. The paper will be constructed as follows: in 
Section ^ we will present the one loop correction to the continuum current evaluated between 
an incoming b quark and outgoing c quark, retaining all terms to order -p-. This is matched to 
continuum Euclidean NRQCD to 0{jj), giving a basis of continuum NRQCD operators for each 



of the currents Vo, 14, Aq, A^. In Section [II we select a basis of operators in the Lattice theory 
and define the one-loop mixing matrix, giving the projection of each element of a basis for a given 
Lattice current at one loop on the basis of continuum NRQCD operators. Finally, we use these 
results to construct relativistic MS currents to this order from lattice operators in Section 



II. CONTINUUM CALCULATION 

We require the continuum renormalisation of the Vector and Axial currents involving an incoming 
b quark with momentum p and outgoing c quark with momentum p' , i.e. Uc{p')Tuiy{p), where either 
r = 7^ or F = 7^1 75- At one loop the graphs in Figures |l] to |^ contribute. The calculation was 
performed in Minkowski space using dimensional regularisation in Feynman gauge, the on-shell MS 
scheme, and having introduced a fictitious gluon mass A to regulate infra-red divergences. The 
treatment allowed for non-zero external momenta, keeping terms linear in and jj . 

In general, we write the perturbative expansion of some quantity, Q, as 



Q 



k = 



We write the Feynman graph in Figure |l| as OsA and the leg corrections in Figures |^ and ^ 
as ^OsTZl^^^ and ^asTZ^^^ respectively. Dropping the external spinors, the total current correction 
may be written as 



AF : 
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Here both A""' ' , and contain (cancelling) infra-red divergences of modulus ^41og -jj. The result 
is IR-finite, and we obtain that the MS current to one loop and 0{jj) is as follows 



where, 
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and 



It is worth noting here that at higher orders in the -^expansion the coefficients d become functions of the recoil lu 
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-1 ;r = 7''75 ' \ 75 ;r = 7"75 

The Gordon decomposition of the spinor gives that, for rric — mi,, 

- M (p + pT - 

wy u = uu 

' 2m 

and so the Vector current between degenerate quarks is, of course, not renormalised. 

In order to match these currents with Euchdean lattice NRQCD currents we first analytically 
continue this resuH to Euchdean space, then apply the Foldy-Wouthuysen expansion of the external 
spinors [B to the appropriate order, taking care with the Euchdeanised 7-matrices. 



A. Analytic Continuation 

We define the Euclideanisation as follows. Denoting the Minkowski gamma matrices and four 
vectors with a caret, 7, so that {7^,71/} ~ ig^v, we use the Dirac-Pauh representation: 
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70 = 
-7j = 
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We also define the Euclidean four vectors in terms of the contravariant and covariant Minkowski 
vectors as: 



X 



,0 



Jb j J 

We introduce Euclidean 7 matrices such that {7^,71^} = 25 
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For consistent notation in Euclidean space, we define Euclidean bilinear currents in terms of the 
Euclidean 7 matrices, so that 

AP = Ao = A° = io 



A, 



iA> 



V = Vj = -iV' = iVj 

Dropping terms in A^ which are subsequently of higher order in , and suppressing the spinors, 
this gives us in Euclidean space: 
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B. Non-Relativistic Expansion 



We now perform the Foldy-Wouthuyseii expansion of the external spinors Uh to order — in terms 
of the Euchdean 7 matrices to produce two component spinors Uq: 
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2m 



1-P 
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where 



UQ 



Uq 





and Uq is a two component spinor, and m is the pole mass of the quark arising from the on- 
shell condition (more formally, m is the one loop pole mass whenever it appears at tree level since 
we are working consistently to order one loop). Inserting this in the results for the continuum 
renormalisation and keeping terms to order results in the following mixed continuum currents, 
with the relevant two component spinors suppressed. 
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We therefore define continuum Euclidean operators 57/ in Table | 
corresponding set of lattice operators, and use the notation 



to which we wish to match a 
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1 — loop 
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The vector current is a conserved current for mt, = mc and, for this reason, the continuum 
renormalisation ( factor in square brackets in ( p^ ) is unity in the ^ — » 1 limit. The operator corre- 
sponding to the temporal component of the vector current, Ipauu, is a conserved current of both the 
lattice NRQCD and HQET actions, and we shall see that the corresponding lattice renormalisation 
is also trivial in the mb = mc limit. 

Unlike the lattice correction, the pure IfQET contribution to the vector current renormalisation 
at zero recoil remains trivial with non-identical quark masses at leading order in — and — . As a 
result the above expressions (p4) and (M) are the same as the renormalisation at zero recoil for the 
continuum heavy quark effective theory, commonly referred to as rjv and rjA respectively Q. This 
is somewhat fortuitous since it will allow us, in section M, t o use the existing result for the BLM 
scale n for the processes contributing to expressions (p^ and with nv ^ 0.92y/mbmc and 
fiA — 0.51y^mbmc. 



III. LATTICE CALCULATION 



A. Lattice action, currents, and Feynman rules 

We use the usual Wilson plaquette action defined in terms of the usual SU(3) gauge links 
f/p(a;), and the tadpole improved O(jg-) lattice NRQCD action [^, 



a£NKQCD = ^HxMx) - ^t(, + at) (1 - ^) (1 - (1 _ i^)" (1 _ ^) 

where 
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and 



Pc^fiix) = -i^ f/c(x)f/^(a; + a)(/J,(x + P)Ulix). 

In the final results quoted in this paper, ug is taken from the plaquette and to one loop order in 
perturbation theory Uq^' ~ ~f • However, sufficient information will presented to allow a reader to 
change these results for different values of uq. To this order in Us any dynamical quarks present in 
the simulation do not appear in the perturbation theory. This calculation is therefore applicable to 
both quenched and dynamical simulations. It is well known how to perform perturbative calculations 
of the renormalisation of weak currents in NRQCD jlo|, and the Feynman rules in the Euclidean 
lattice theory used here are the same as those published by Morningstar and Shigemitsu This 
paper completes a programme of work previously discussed at two lattice conferences [ pT| . 

Inspection of the continuum currents suggests the basis for the lattice current operators given in 
Table The momentum dependent operators are reproduced on the lattice with gauge covariant 
derivatives, 

VMX(a;) = ^ [U^{x)x{x + A) - Ul{x - fi)x(x - fl)] 

which introduces additional quark gluon vertices that are not present in the continuum theory, as 
well as a tadpole improvement counter term. For convenience in simulations, we define the lattice 
operators in terms of the bare mass, Mq. Since we require that the one-loop renormalised kinetic 
meiss in NRQCD matches the one-loop renormalised mass in the relativistic theory |^ we must 
introduce additional Zm factors for lattice currents wherever the bare mass Mq appears at tree 
level. The mass ratio ^ appears in the continuum coefficients only at order a^, and hence the bare 
ratio ^ can be used to this order. 

We have additional Feynman rules associated with the current operators involving derivatives 
(and hence gauge links), giving both one and two gluon couplings. The conventions for momenta 
involved in current insertions of an operator O both at leading order, and couplings to both 1 and 
2 gluons of polarisation j are as given in figure |^. 

The Feynman rules associated with the spatial vector and temporal axial currents at the vertex 
are simply the insertion of the appropriate Pauli structure. The Feynman rules associated with the 
spatial vector current insertions are as follows 



O]^" (tree) = - ^ sin fcjcrfccr. 
Ol" (1 - gluon) = - ^ cos [k, - ig, 
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<M2 - gluon) = ^ sin 



akfJj 



z invariance implies that the kinetic and pole masses are consistent in the relativistic theory 



03*= (tree) 
C)J''=(l-gluon) 
03^M2-gluon) 

where the 1 and 2 gluon couphngs are to gluons of polarisation j. 

The rules associated with the leftward derivative currents are both the obvious analogues of the 
above and not in fact required since the diagrams involving the second current may be switched 
for those of the first current with the masses swapped and momenta reversed, halving the required 
coding effort. The Feynman rules associated with the first temporal axial lattice current are as 
follows, 

of « (tree) = ^ ^ sin k,a, 
j 

Ofil - gluon) = ^ cos (fc, - ig,) aj 
of" {2 - gluon) = ^ sin (fc, - i(gi, + qa,)) a, 
while, as above, the rule for the second current is not required. 
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B. Mixing Matrix 

We define the 1-loop mixing matrix of the lattice theory for each of the four currents J G 
{Vb, Vk,Ao, Ak} as follows: 



where the range of the summation over (3 depends on the dimension of the basis for each continuum 
current. To avoid ambiguity, we shall denote the symbols for lattice Feynman diagrams in bold 
font to differentiate them from symbols used in the continuum calculation (although the different 
symbols do not in fact appear in the same equations). For both the temporal vector and spatial 
axial currents the mixing matrix is rank one, and we have (where J is either Vo or Ak): 



zj^' gives the relevant leg correction, whose relation to the self energy diagram (Figure P) is described 
in appendix Q where we have followed closely the conventions of Morningstar . A'''^' corresponds 
to contributions from the Feynman diagram of Figure ^ where gluon couplings from each of the Vt, 
^— , and terms in the action are included at each vertex. 

For the spatial vector current (and indeed the temporal axial current) the mixing matrix is non- 
trivial. We denote the vertex correction to each of the lattice currents O^'' from the sum of Feynman 
diagrams represented in Figures Hand I by aI*^'^'. At one loop the the vertex correction overlaps 
with each of the corresponding continuum operators, giving rise to the non-trivial mixing problem 

We resolve the Pauli and momentum structure of aJ^*"'^' using the following projections with j 7^ k, 
then setting the spatial momenta to zero: 

Ar^'^l = 2iA/,a,a.4Ar^-'y Ar^'^' = 2iM.aka, ^A^^^'^ 

iVj3 — inib-7^J\.-^ 2-'^ll ■''■14 — ^^'^cg^^Vj^ 2-'^12 ■ 



We may use the bare lattice mass for convenience in place of the the pole mass in these projections 
since we are evaluating one-loop coefficients, the difference being 0{as)- Mixing matrix elements 



A.2j3 are related to Ai^ by Mj, ^ Mc, and it will be seen later that Aa^ and A4^ are not required 
for our calculation. Hence we write the contributions to the mixing matrix for the spatial vector 
current as follows, 



•^r?3 = ^'='13 + 
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where the the Tim factors arise from matching the bare mass in the lattice current to the pole mass 
in the continuum current to order Us, and the tadpole improvement counter terms are 



yTI 

.Tl^k [1 
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Similarly, for the temporal axial current we project the vertex corrections on to the continuum basis 
using, 



There is also a tadpole improvement counter term for the temporal axial current, 

.Ti^ori] „ [1] 

We write the contributions to the mixing matrix for the temporal axial current as follows. 



M 



Q/3 
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+{ZW + Z^i'i}5,25./3 + A-tf^ + A 



/3 

TI-*o [1] 
■-a/3 



The major task in the renormalisation is to calculate the terms in the above mixing matrices nu- 
merically using the Feynman rules for the lattice action. 

The integrands for each of the Feynman graphs (and their first and second derivatives with respect 
to both Pi and p'j where appropriate) were coded in two independent ways in order to verify the 
correctness of the calculation. 

A C-|— I- class was written to encode the Pauli algebra, and wrapped in a Taylor series class giving 
the first few entries of the (double) expansion in pj and p^. The standard operators were overloaded 
and standard functions applied the chain rule to their arguments. This enabled the Feynman 
diagrams to be constructed easily from the Feynman rules, with the Pauli manipulations taken 
care of by the code. The leading entries in the double Taylor series gave the first few derivatives 
with respect to Pj and p'j.. Additionally much faster Fortran codes were written with the Pauli 
manipulations carried out manually, and all derivatives with respect to the external momenta taken 
by hand. VEGAS was used to evaluate the one loop integrals. 

The lattice effective theory is constructed to have the same infra-red physics as the relativistic 
continuum theory, so that we find the same (cancelling) logarithmic infra-red divergences in the 
lattice renormalisation in the lattice diagrams as were present in continuum theory. These arise 
from graphs where a temporally polarised gluon is exchanged via the Vt coupling in the action. 
To handle the infra-red divergent graphs we always perform the integral over fco analytically, and 
then integrate infra-red finite sums of graphs, e.g. ^{z|^' -I- z|^|} -I- A^\', over the remaining three 
dimensions corresponding to the spatial components of the loop momentum. 



IV. MATCHING AND RESULTS 



We wish to form a linear combination of the lattice currents such that they match the continuum 
currents at one loop. We reconstruct the MS continuum currents from the lattice currents as follows, 

= f: [Zr^+aspj] {c{p')\Oi\h{p)Y-'°°^ 
where Z'-j^"" G {0, 1} is the tree level piece of Zf^ 
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In what follows we shall present results for the for each of the currents in turn as a function 
of both the bare lattice mass and the mass ratio £, = ^ . For reasons of practicality, we make 
the stabilisation parameters for each quark, and nt, implicit functions of the mass as indicated 
in table 



III 



A. Temporal Vector Current 



For the temporal vector current we have 





= {1^ 
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^ 3^ 


-3i±|logC-6 





We plot the renormalisation coefficient in figure In table ^ we give the correction as a 
function of -jj- for ^ = 0.1,0.2,0.3,0.4,0.5 and 1.0, useful ratios for lattice simulations. Since the 
lattice current is conserved, the one-loop renormalisation coefficient vanishes, for all Af;, in the case 
^ = 1, as clearly shown by figure 

Luke's theorem shows that there is no 0(-p) correction to form factors involving the vector 
current in HQET at zero recoil. Luke's theorem is a non-perturbative observation based upon the 
symmetries of HQET. The required symmetries are satisfied by Lattice NRQCD, and the theorem 
should certainly be manifested in the one-loop coefficients at zero-recoil. Demonstrating this is less 
trivial than in HQET however, since M dependence arises out of terms in the action as well as 
through current insertions, and we proceed by numerically calculating the M dependence of the 
lattice NRQCD contribution to the renormalisation. In figure ^ we display the lattice one loop 
renormalisation — plotted versus (-p for various values of ^ = This plot shows 



that the leading correction is linear in ( 



Me 
1 
Mc 



Mi, 



in the 



_i_ 

M 



region, with a universal slope. 



This is consistent with the absence of the -jg- term, as expected by Luke's theorem applied to lattice 
NRQCD. 

Of course the analogue of Luke's theorem should be held non-perturbatively by lattice NRQCD, 
and a fit of the dependence of the non-perturbative lattice form factor on -jg- near the static limit 
should be made on simulation data to check the consistency with Luke's theorem. 



For the spatial axial current. 



B. Spatial Axial Current 





= {l + a,p■*'=}o■^^ 
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We plot the renormalisation coefficient p^*' in figure ^ and give the correction as a function of 
for I = 0.1, 0.2, 0.3, 0.4, 0.5 and 1.0 in table 0. 
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C. Temporal Axial Current 

For the temporal axial current we have 

= { 1 + « } Of « + { 1 + a.pf « } Of 



where 
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and 



-3i±|logC-8-6i±|-12^1°sC 



loge-8 + 6ii| + 12 



We plot the renormalisation coefficients p, ° and ° in figures |l| and Tile various contributions 
to pf° for C = 0.3 are given in table Knl allowing the reader to see relatively easily the change 
that must be applied to obtain the coefficients with a different tadpole improvement prescription, 
such as the mean link in Landau gauge. The diagonal elements App of the mixing matrix both 
contain infra-red divergences which exactly cancel those of the wavefunction renormalisation, and 
have associated tadpole improvement counter terms associated with the uo factor in the derivative 
operator. It is, in fact, interesting to observe the tadpole improvement programme in action; there 

TIfll fll TI'^'^ [ll 1 fll 

are large cancellations between the Z™ and Z™ terms and between the A^^ and the f {Z|^ + 
terms. This cancellation is even more manifest when we isolate the tadpole graph 



Ao[l] 



contributions within Afj"'^' and Zm. The two overall co rrect ion coeffi cients are tabulated as a 



function of -j^ ior ^ = 0.1, 0.2, 0.3, 0.4, 0.5 and 1.0 in tables |VIl| and |VII] 



D. Spatial Vector Current 

For the spatial vector current we obtain 

VF^ = {1 + a.pr" } + { 1 + ^sP^" } O^" + a.pl^ Ol" + a.pl^ 0\- 

where 
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Since the O3 and ©4 do not appear at tree level, we can see that it is in fact unnecessary to calculate 



and A 



. p2' in figures | 



.5 to h 8l The various 



1.3^ diiu J 14^. We plot the renormalisation coefficients pj 
contributions to p^^ and p^*" are given for ^ = 0.3 in table ^ and table The four overall 
co rrecti on coefficients are given as a function of ^ for ^ = 0.1, 0.2, 0.3, 0.4, 0.5 an 1.0 in tables |Xl| 

to Efvl. 



V. RELEVANT MOMENTUM SCALE FOR as 



In section ^ we argued that the regularisation scheme correction is dominated by ultra-violet 
physics and is therefore calculable within perturbation theory. This is addressed in this section by 
attempting to calculate the appropriate scale for as in the correction to Vb and Ak- 

The Lepage-Mackenzie q* scale ^ for some one-loop process whose Feynman amplitude A = 
I (l^-'^(^) defined by 



iog(agT = - "T (77) 



/ ^J(fc)log(a fc)^ 



When a one- loop expression is made up of two sub-processes, A = J ^^^^4 lA^k) and B 
I (2ir)'^ *® nis-y combine the q*'s as follows, 



/ (iF^^(fc) l«gK)' + / (l^^s(fc) log(afc) 

iq — 

or more conveniently. 



iog(ag )a+b = f rf^fe r f m ^^^^ 



ql+B = [{gl)^ X {qif]^ . (79) 

Here we shall calculate the Lepage-Mackenzie q* from lattice per turb ation theory, and then use this 
rule to combine it with the continuum BLM scales from section [[IB| to obtain an overall scale. 

The combination of lattice graphs containing the cancelling IR divergences proved troublesome 
using 4 dimensional numerical integration, so we used the (non-rigorous) Hernandez and Hill pre- 
scription to deal with the kg integral , an d integrated the remaining spatial momenta numerically. 



We present the results for qy^ in table XV . The overall scale is reassuringly large, and is dominated 
by the continuum scale due to the typical relative size of the corrections. We therefore believe that 
the use of Hernandez and Hill prescription to obtain numerical stability constitutes an irrelevant 
error in the overall scale. We note that av{q*) < 0.22 for q* > 3GeV, which is the case in a typical 
simulation regimes for Vb |lg|. 



We present the results for in table XVl . There is a line (straddled by the bold highlighted 
data in tab le XVl ) in the ^ — Mb plane where q* becomes poorly behaved. It can be seen in 
table XVII that the individual scales g^'^* and Qa™"* Eire large at this point, but that the overall 
one-loop correction is near vanishing. The location of this uncontrolled behaviour is not dependent 
on our use of the Hernandez and Hill prescription, however the details of the behaviour will, of 
course, be somewhat dependent. The zero in the overall correction introduces a vertical asymptote 
for log(a(j'*)^, due to the vanishing denominator in ( |7^ ) and consequently a limit of infinity or zero 
for q* depending on the direction of approach. The resultant low scale on one side of the zero in the 
correction is not considered problematic by the authors, since it is a breakdown of the mean scale 
interpretation of ([77[), that occurs whenever the correction is smaUj^ 

The prescription used to set the scale at which to evaluate as in a mixing calculation is considered 
an open question by the authors, however we remain encouraged by the reliability of the perturbation 
theory suggested by the scales determined thus far. 



VI. SUMMARY 

We have shown how to construct the MS weak decay currents for B ~* D'*' transitions from 
operators in the lattice NRQCD effective theory at and around the zero recoil point for a wide range 
of p henomenologically relevant lattice mass parameters. These constants are of use to both recent 
|lq] and future [hsl calculations of the semileptonic decay form factors. The one-loop coefficients 
typically have modulus < 0.3 in the usual simulation regime. The lattice temporal vector current 



*As noted by Morningstar the integrand lA{k) + Isik) must be either strictly positive or strictly negative for the mean 
value theorem to guarantee < aq* < 2ti. Exact cancellation of the integrals is likely to exacerbate the problem. 
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is not renormalised in the degenerate quark case (this is not the case for Wilson fermions), with 
the result that the current correction for physical ^ ratios is particularly small. For two of the 
currents we have estimated the Lepage-Mackenzie q* scales, and find that it is large, giving overall 
one-loop corrections of order < 5%, and suggesting perturbation theory is converging adequately. 
With further suppression, the next (and higher) order corrections in the operator renormalisation 
could well be below the expected statistical error of the BaBar data. |^ 



VII. ACKNOWLEDGMENTS 

We wish to thank Junko Shigemitsu for many useful conversations, and cross checking the mass 
renormalisation. We made use of unpublished notes by Colin Morningstar for checking our contin- 
uum calculation. Peter Boyle is supported by PPARC grant PP/CBA/62 . Both Peter Boyle and 
Christine Davies fondly acknowledge the hospitality of the UCSB, Santa Barbara where some of this 
work was carried out. 



[1] The BaBar Physics Book: Physics at an asymmetric B factory. 

By BABAR Collaboration (P.P. Harrison, e d. et al.). SLAC- R-0504. Oct 1998. 



[2] A. Kronfeld, FERMILAB-CO NF-98-393-T, |hep-ph/9812288 



L. LeUouch, hep-ph/9912353 
S. Aoki, tiep'ph/9912lif 



S. Hashimoto, Nucl. Phys. (Proc Suppl), 83-84, 3, (2000). 
[3] L. Foldy and S. Wouthuysen Phys. Rev. D,48, 29 (1950) 
[4] M. Neubert Phys. Kept, 245, 259, (1994). 

[5] S. Brodsky, CP. Lepage, P. Mackenzie Phys. Rev. D 28, 228, (1983). 
[6] M. Neubert Phys. Lett. B, 341, 367, (1995). 
[7] K. Wilson Phys. Rev. D, 10, 2445 (1974). 

[8] G.P. Lepage and P. Mackenzie Phys. Rev. D 48, 2250, (1993). 
[9] B. Thacker and P. Lepage, Phys. Rev. D 43, 196, (1991). 

P. Lepage et. al. Phys. Rev. D 46, 4052, (1992) . 
[10] C. Morningstar, J. Shigemitsu 

Phys. Rev. D, 57, 6741, (1998). 

Phys. Rev. D, 59, 094504, (1999). 

Nucl. Phys. (Proc Suppl), 63, 341, (1998). 
[11] P. Boyle and C. Davies, Nucl. Phys. (Proc Suppl), 73, 378, (1999). 

P. Boyle and C. Davies Presented at Lattice 99, Pisa, July 1999 
[12] C. Morningstar, Phys. Rev. D 48, 2265, (1993). 
[13] G.P. Lepage J. Comput. Phys. 27 192 (1978). 
[14] M.E. Luke Phys. Lett.B 252, 447 (1990). 
[15] O. Hernandez and B. Hill, Phys. Rev. D,50, 495, (1994) 



[16] J. Hein Presented at Lattice 99, Pisa, July 1999. |hep-lat/990805 j 
[17] C. Morningstar, Phys. Rev. D 50, 5902, (1994). 
[18] C. Davies, J. Hein, et al tn preparation. 

[19] A. Czarnecki, K. Melnikov Phys. Rev. Lett. 78, 3630, (1997). 



APPENDIX A: SELF ENERGY CALCULATION 

To one loop the graph topologies shown in figures ^ and ^ contribute to the lattice self energy. The 
diagrams contain contributions from both the temporally polarised gluon vertex, and from spatially 



Certainly, the two loop renormalisation correction to HQET operators is known to be 0(1%) [M. 
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polarised gluons via the ^ term and via the a ■ B term. The lattice self energy must satisfy cubic 
invariance in the spatial momenta, and we write the first terms of a double taylor series for the 
self-energy as 

T,^'"\apo,ap) = g"^ \ flo + ^liapo + O.2 Jl " + ■■■> (Al) 
I lAloa J 

where clearly, 

/!^o = S(0,0) (A2) 

= -i-^S(0,0) (A3) 
aapo 



=aMo^^E(0,0) (A4) 

(A5) 

2 

to order a then the small ^ expansion of the heavy quark propagator is 

G-\p) = G-\p) - T.{apo,ap) (A6) 

, Po I ■ 2o • 2o 2o I KT, 

2 2 —2 2 

2 2 

where to order g and ^ we have, 

jpo = ipo - g^^o (A9) 
^^ = l-/(!^o + ni) (AlO) 
^5;^" = l + ff'(t^2-ni) (All) 

The contributions from the tadpole diagram figure ^ to Q,o and Sli explicitly cancel in due to 
the equation of motion of the external spinor. In the Zm calculation, both the tadpole and normal 
topologies contribute. In order to compute Z,n it was necessary to code a routine to compute 
the second derivative of the self energy with respect to one of the external spatial momentum 
components. This was fairly involved, due to the nature of the lattice Feynman rules. There are 
cancelling infra-red divergences in both Qq and in Qi , so that the integral over the ko was performed 
analytically and the remaining three dimensions integrated numerically. Our numbers for both the 
Z^, Ae, and Zm agree within statistical error with those of Morningstar and Shigemitsu px| ]. 



1 d' 



,2 
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p[l] 

FIG. 1. QfaA : Continuum vertex correction to & — + c scattering diagram via weak current F = 7;j or F = 7^75. 



r 




Uc(p',mc) Ub(p,mb) 





FIG. 4. Feynman rules associated with insertion of a lattice weak current operator to order g", , and g^. The gluon 
couplings only arise for T = Vi and T = Ao- Here (3 runs over the dimension of the lattice operator basis corresponding to the 
continuum current F, as given in Table 
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FIG. 5. Regular topology self energy contribution. 
Contributions arise from a temporally polarised gluon coupling to vertices arising from the dt term in the action, and spatially 
polarised gluons coupling to vertices arising from the ^ term and the term. For the wavefunction renormalisation the 

^'^'^ 2^ coupling diagram vanishes at p = 0, but, since the derivatives do not, the mixed diagrams contribute to the 
kinetic mass renormalisation. 

k 




P P+k P 
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FIG. 6. Regular topology vertex correction. 

2 

Temporally polarised gluons couple through the dt term in the action for and spatially polarised gluons through both the ^ 
and term. Graphs involving mixed ^ and couplings vanish for the temporal vector and spatial axial currents. 




FIG. 7. Diagrams of this topology contribute to both the temporal axial and spatial vector currents. The "tadpole" arises 
from gauge link in the derivative in the current operators. 

k 




FIG. 8. Diagrams of "left ear" and "right ear" topologies contribute to both the temporal axial and spatial vector currents. 
The gauge link in the spatial derivative in the current operators creates a coupling to spatially polarised gluons, and the 

2 D 

coupling to the leg arises from both the ^ and terms. 
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FIG. 9. Tadpole contribution to the self energy. 
Due to the external equation of motion, this diagram does not in fact contribute to the wavefunction renormalisation, however 

2 o 

it does contribute to the kinetic mass renormalisation, through both ^ and two gluon couplings. 



i f 



FIG. 10. p^o as a function of ^ and jjl^. 




FIG. 11. Pure lattice correction — as a function of — -jg^)^ for curves ^ = 0.2 . . . 0.5 running from top to bottom. 
The leading correction is linear and universal showing that Luke's theorem is held to one-loop in lattice NRQCD. 
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FIG. 14. as a function of £ and — . 
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FIG. 16. p^*" as a function of ^ and 



0.5 r 




FIG. 17. p^'' as a function of ^ and 
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FIG. 18. pX'' as a function of (, and — . 




TABLE I. Continuum operators for which we will select corresponding Lattice operators. 











^l"^ — Ipauli 


nf = —iau 


n*''= — —irr, "'^ 
4 rric 


CtAo _ a-p 
"l - 2m. 
O'*o _ CT-p 
"2 - 2mc 
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TABLE II. Bases of Lattice Currents to be matched to each of the weak currents. 
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^\ ^ — Ipauli 


Ot" = -iak 




— 2M<: " 
_ Vfc 

_ Vfc 

— Me 


of" 


- -i^ 



TABLE III. Selected stabilisation parameter in different mass regions. 



M > 0.8 


M > 1.0 


M > 1.2 


M > 1.6 


M > 4.0 


n = 5 


n=4 


n=3 


n=2 


n=l 



TABLE IV. The one loop piece of the temporal vector current renormahsation for ^ = 0.1, 0.2, 0.3, 0.4, 0.5, and 1.0 . 



rrib 






P-*MC = 0.3) 


P^HC = 0.4) 


P-^M^ = 0.5) 


P^H^ = m 


2.000 








0.127146(7) 


0.063266(3) 





3.000 






0.168931(8) 


0.084085(3) 


0.043936(2) 





4.000 




0.28746(1) 


0.131438(5) 


0.068204(2) 


0.0366150(9) 





6.000 




0.208433(6) 


0.103350(2) 


0.0564405(10) 


0.0311405(5) 





8.000 


0.46550(2) 


0.178367(4) 


0.092615(1) 


0.0516583(6) 


0.0288833(3) 





10.000 


0.39874(1) 


0.163900(3) 


0.0871368(10) 


0.0492116(4) 


0.0276885(2) 





12.000 


0.360656(8) 


0.155447(2) 


0.0839188(7) 


0.0477178(3) 


0.0270377(1) 





16.000 


0.320742(5) 


0.146303(1) 


0.0803731(4) 


0.0462188(2) 


0.02633675(8) 






TABLE V. The one loop contribution to the spatial axial current renormahsation for ^ = 0.1, 0.2, 0.3, 0.4, 0.5 and 1.0 . 





P^H? = o.i) 


P^H? = 0.2) 


P^H$ = 0.3) 


p^'=(^ = 0.4) 


p^M^ = 0.5) 


p^''(^ = 1.0) 


2.000 








0.23875(1) 


0.117205(8) 


-0.066846(4) 


3.000 






0.161977(8) 


0.031141(5) 


-0.036802(4) 


-0.140394(2) 


4.000 




0.251373(9) 


0.042554(5) 


-0.048200(3) 


-0.096478(2) 


-0.169133(1) 


6.000 




0.081220(4) 


-0.049163(2) 


-0.109181(2) 


-0.142730(1) 


-0.1917372(5) 


8.000 


0.346311(8) 


0.016752(3) 


-0.083899(1) 


-0.1327985(9) 


-0.1605718(7) 


-0.2001877(3) 


10.000 


0.248523(5) 


-0.014434(2) 


-0.1012549(9) 


-0.1445488(6) 


-0.1695061(4) 


-0.2043054(2) 


12.000 


0.192853(4) 


-0.032424(1) 


-0.1112436(7) 


-0.1513905(4) 


-0.1744741(3) 


-0.2066181(1) 


16.000 


0.134707(2) 


-0.0515471(7) 


-0.1218945(4) 


-0.1583394(3) 


-0.1796509(2) 


-0.20898845(8) 



TABLE VI. Contributions to the temporal axial current renormahsation for ^ = 0.3 . 



rrib 


rib 


rric 


Tic 


i{Z^,H-Z^J + Afi« 




Bi 




7TI 


\Ao TI 


Pi 


3.0 


2 


0.9 


5 


-1.3382 


-0.6789 


-0.3809 


1.1458 


-0.7854 


1.0472 


0.2286 


4.0 


2 


1.2 


3 


-1.2009 


-0.5833 


-0.3809 


1.0861 


-0.8508 


1.0472 


0.1209 


6.0 


1 


1.8 


2 


-1.0786 


-0.4599 


-0.3809 


0.8770 


-0.7854 


1.0472 


0.0189 


8.0 


1 


2.4 


2 


-1.0188 


-0.3635 


-0.3809 


0.8761 


-0.8508 


1.0472 


-0.0710 


10.0 


1 


3.0 


2 


-0.9846 


-0.3017 


-0.3809 


0.8749 


-0.8901 


1.0472 


-0.1266 


12.0 


1 


3.6 


2 


-0.9623 


-0.2574 


-0.3809 


0.8733 


-0.9163 


1.0472 


-0.1654 


16.0 


1 


4.8 


1 


-0.9355 


-0.1990 


-0.3809 


0.8708 


-0.9490 


1.0472 


-0.2153 


32.0 


1 


9.6 


1 


-0.8989 


-0.1040 


-0.3809 


0.8673 


-0.9981 


1.0472 


-0.2944 


64.0 


1 


19.2 


1 


-0.8806 


-0.0534 


-0.3809 


0.8649 


-1.0227 


1.0472 


-0.3363 


128.0 


1 


:!8. 1 


1 


-0.8721 


-0.0270 


-0.:!809 


0.8():!() 


-1.0:! 19 


1.0172 


-0.:!574 
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TABLE VII. pf° contribution to the temporal axial current ronormalisation. 



rrib 


pf°($ = 0.1) 


pf«(^ = 0.2) 


pf«(^ = 0.3) 


pf«(C = 0.4) 


^^(^ = 0.5) 


pf»(C = i.o) 


2.000 


- 


- 


- 


0.195(3) 


0.070(3) 


-0.138(3) 


3.000 


- 


- 


0.229(3) 


0.068(3) 


-0.025(2) 


-0.159(2) 


4.000 




0.385(3) 


0.121(2) 


-0.009(2) 


-0.082(2) 


-0.169(2) 


6.000 




0.224(2) 


0.019(2) 


-0.075(2) 


-0.126(2) 


-0.159(2) 


8.000 


0.584(2) 


0.094(2) 


-0.071(2) 


-0.141(2) 


-0.173(2) 


-0.172(1) 


10.000 


0.435(2) 


0.008(2) 


-0.127(1) 


-0.181(1) 


-0.203(1) 


-0.180(1) 


12.000 


0.326(2) 


-0.052(1) 


-0.165(1) 


-0.208(1) 


-0.223(1) 


-0.185(1) 


16.000 


0.175(2) 


-0.130(1) 


-0.215(1) 


-0.243(1) 


-0.247(1) 


-0.190(1) 


TABLE VIII. contribution to the temporal axial current renormalisation. 


mb 


P^(^ = 0.1) 


P^(^ = 0.2) 


P^(^ = 0.3) 


P^(^ = 0.4) 


p^(e = 0.5) 


P^(^ = l-0) 


2.000 








-0.027(5) 


-0.112(4) 


-0.142(3) 


3.000 






-0.041(5) 


-0.129(4) 


-0.149(3) 


-0.158(2) 


4.000 




0.067(5) 


-0.095(4) 


-0.134(3) 


-0.145(3) 


-0.168(2) 


6.000 




-0.020(4) 


-0.079(3) 


-0.106(3) 


-0.118(2) 


-0.158(2) 


8.000 


0.239(5) 


-0.005(3) 


-0.052(3) 


-0.077(2) 


-0.104(2) 


-0.170(2) 


10.000 


0.182(4) 


0.020(3) 


-0.025(2) 


-0.063(2) 


-0.070(2) 


-0.179(1) 


12.000 


0.162(4) 


0.043(3) 


-0.006(2) 


-0.027(2) 


-0.066(2) 


-0.183(1) 


16.000 


0.182(3) 


0.083(2) 


0.038(2) 


-0.016(2) 


-0.060(2) 


-0.192(1) 



TABLE IX. Various contributions to the pj*" spatial vector current renormalisation for f = 0.3 . 



nib 


rib 


rric 


Tic 




■''■21 


Bi 




7T1 
"»/) 


A Vfc TI 

■''^11 


Pi 


3.0 


2 


0.9 


5 


-1.4361 


0.6402 


0.2873 


1.1496 


-0.7854 


1.0472 


-0.3281 


4.0 


2 


1.2 


3 


-1.2605 


0.5574 


0.2873 


1.0932 


-0.8508 


1.0472 


-0.2991 


6.0 


1 


1.8 


2 


-1.1096 


0.4418 


0.2873 


0.8795 


-0.7854 


1.0472 


-0.1862 


8.0 


1 


2.4 


2 


-1.0372 


0.3520 


0.2873 


0.8786 


-0.8508 


1.0472 


-0.1024 


10.0 


1 


3.0 


2 


-0.9971 


0.2933 


0.2873 


0.8748 


-0.8901 


1.0472 


-0.0407 


12.0 


1 


3.6 


2 


-0.9719 


0.2510 


0.2873 


0.8734 


-0.9163 


1.0472 


0.0039 


16.0 


1 


4.8 


1 


-0.9406 


0.1951 


0.2873 


0.8720 


-0.9490 


1.0472 


0.0626 


32.0 


1 


9.6 


1 


-0.8983 


0.1031 


0.2873 


0.8696 


-0.9981 


1.0472 


0.1638 


64.0 


1 


19.2 


1 


-0.8806 


0.0527 


0.2873 


0.8663 


-1.0227 


1.0472 


0.2243 


128.0 


1 


38.4 


1 


-0.8721 


0.0268 


0.2873 


0.8649 


-1.0349 


1.0472 


0.2555 



TABLE X. Various contributions to the pg*" spatial vector current renormalisation for ^ = 0.3 . 



nib 


rib 


rric 


ric 






B3 


p-.i 


3.0 


2 


0.9 


5 


-0.0091 


-0.6212 


-0.1467 


0.4836 


4.0 


2 


1.2 


3 


-0.0052 


-0.5416 


-0.1467 


0.4001 


6.0 


1 


1.8 


2 


-0.0030 


-0.4333 


-0.1467 


0.2896 


8.0 


1 


2.4 


2 


-0.0017 


-0.3472 


-0.1467 


0.2022 


10.0 


1 


3.0 


2 


-0.0010 


-0.2897 


-0.1467 


0.1439 


12.0 


1 


3.6 


2 


-0.0007 


-0.2483 


-0.1467 


0.1022 


16.0 


1 


4.8 


1 


-0.0004 


-0.1933 


-0.1467 


0.0470 


32.0 


1 


9.6 


1 


-0.0001 


-0.1024 


-0.1467 


-0.0442 


64.0 


1 


19.2 


1 


-0.0000 


-0.0527 


-0.1467 


-0.0940 


128.0 


1 


38.4 


1 


-0.0000 


-0.0267 


-0.1467 


-0.1200 
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TABLE XI. p^'' contribution to spatial vector current renormalisation. 



mt 


p^''{i = 0.1) 


pr'=(C = 0.2) 


prMe=o.3) 


Pi'i^^OA) 


prMc=o.5) 


Pi'iC^l-O) 


2.000 


- 


- 


- 


-0.364(3) 


-0.396(3) 


-0.339(3) 


3.000 


- 


- 


-0.328(3) 


-0.346(2) 


-0.329(2) 


-0.225(2) 


4.000 




-0.247(3) 


-0.299(2) 


-0.279(2) 


-0.250(2) 


-0.149(2) 


6.000 




-0.214(2) 


-0.186(2) 


-0.151(2) 


-0.121(2) 


-0.034(2) 


8.000 


-0.102(2) 


-0.142(2) 


-0.102(2) 


-0.073(2) 


-0.050(2) 


0.018(2) 


10.000 


-0.095(2) 


-0.077(2) 


-0.041(2) 


-0.016(1) 


0.000(1) 


0.055(1) 


12.000 


-0.062(2) 


-0.025(2) 


0.004(1) 


0.022(1) 


0.035(1) 


0.079(1) 


16.000 


0.009(2) 


0.047(1) 


0.063(1) 


0.072(1) 


0.079(1) 


0.110(1) 


TABLE XIL p^'' contribution to spatial vector current renormalisation. 


nil, 


/4''(4 = ().i) 


(C = 0.2) 


pV' = (>■■'') 


/4'(^ = o.i) 


/4"(i; =()..-,) 


P,"{^ = 1.0) 


2.000 








-0.168(5) 


-0.266(4) 


-0.328(3) 


3.000 






-0.148(4) 


-0.240(4) 


-0.250(3) 


-0.216(2) 


4.000 




-0.026(5) 


-0.181(4) 


-0.208(3) 


-0.200(3) 


-0.139(2) 


6.000 




-0.089(4) 


-0.119(3) 


-0.113(2) 


-0.097(2) 


-0.029(2) 


8.000 


0.145(5) 


-0.056(3) 


-0.055(2) 


-0.040(2) 


-0.031(2) 


0.026(1) 


10.000 


0.090(4) 


-0.009(3) 


-0.002(2) 


0.006(2) 


0.031(2) 


0.060(1) 


12.000 


0.079(4) 


0.029(2) 


0.037(2) 


0.058(2) 


0.063(2) 


0.083(1) 


1(). ()()() 


().i()<)(:5) 


().()<);i(2) 


().i()7(2) 


l).l().")(2) 


0.101(2) 


l).ili(l) 


TABLE XIII. p}^'' contribution to spatial vector current renormalisation. 


nib 






P3^'=(^ = 0.3) 


prn^=o.4) 


P3^''(^ = 0.5) 


P3''''(^=l-0) 


2.000 








0.487(2) 


0.418(2) 


0.2193(9) 


3.000 






0.484(2) 


0.386(1) 


0.315(1) 


0.1435(7) 


4.000 




0.545(2) 


0.400(1) 


0.305(1) 


0.2393(9) 


0.0955(5) 


6.000 




0.433(1) 


0.290(1) 


0.2045(8) 


0.1519(7) 


0.0441(4) 


8.000 


0.592(2) 


0.329(1) 


0.2022(8) 


0.1340(7) 


0.0927(6) 


0.0105(3) 


10.000 


0.510(1) 


0.2526(9) 


0.1439(7) 


0.0875(6) 


0.0529(5) 


-0.0109(2) 


12.000 


0.440(1) 


0.1970(8) 


0.1022(6) 


0.0540(5) 


0.0255(4) 


-0.0256(2) 


16.000 


0.329(1) 


0.1214(6) 


0.0470(4) 


0.0104(4) 


-0.0102(3) 


-0.0446(2) 


TABLE XIV. pp° contribution to spatial vector current renormalisation. 


rub 


pr''(e = 0.1) 


P?(e = 0.2) 


P?(^ = 0.3) 


P?-(C = 0.4) 


P?(C = 0.5) 


P?(C = l-0) 


2.000 








0.1458(7) 


0.1664(8) 


0.2188(9) 


3.000 






0.0997(5) 


0.1171(6) 


0.1244(6) 


0.1429(7) 


4.000 




0.0661(4) 


0.0823(4) 


0.0930(5) 


0.0952(5) 


0.0946(5) 


6.000 




0.0478(3) 


0.0532(3) 


0.0517(3) 


0.0501(4) 


0.0438(4) 


8.000 


0.0259(2) 


0.0339(3) 


0.0310(3) 


0.0271(3) 


0.0234(3) 


0.0100(3) 


10.000 


0.0202(2) 


0.0212(2) 


0.0159(2) 


0.0106(2) 


0.0089(2) 


-0.0112(2) 


12.000 


0.0157(2) 


0.0116(2) 


0.0048(2) 


0.0014(2) 


-0.0047(2) 


-0.0258(2) 


16.000 


0.0081(1) 


-0.0015(1) 


-0.0080(2) 


-0.0158(2) 


-0.0224(2) 


-0.0448(2) 
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TABLE XV. The overall Lepage-Mackenzie scale aqy^ for the current Vb as a function of A4t and ^. 



nib 


agvn(e = 0.1) 


ag?>je = 0.2) 


aq?>„(^ = 0.3) 


ag?>n(C = 0.4) 


aq*y„{^ = 0.5) 


ag?>n(C = 1-0) 


2.000 








1.6 


1.6 


oo 


3.000 


- 


- 


1.7 


1.7 


1.8 


oo 


4.000 




1.7 


1.8 


2.0 


2.2 


oo 


6.000 




2.1 


2.4 


2.8 


3.1 


oo 


8.000 


2.0 


2.6 


3.2 


3.8 


4.3 


oo 


10.000 


2.3 


3.3 


4.1 


4.9 


5.5 


oo 


12.000 


2.7 


4.0 


5.1 


6.0 


6.8 


oo 


16.000 


3.6 


5.5 


7.0 


8.3 


9.5 


oo 



TABLE XVL The overall Lepage-Mackenzie scale oQa^ for the current Ak as a function of Mi, and 5. 



nib 


agije = o.i) 




agi, = 0.3) 


aql, (C = 0.4) 


agi, (5 = 0.5) 


agA.(C = l-0) 


2.000 








3.6 


6.8 


0.4 


3.000 






3.1 


27.9 


0.2 


1.6 


4.000 




2.1 


5.1 


0.8 


1.4 


2.2 


6.000 




1.8 


2.1 


2.3 


2.5 


3.4 


8.000 


1.6 


0.5 


3.1 


3.1 


3.3 


4.4 


10.000 


1.6 


24.0 


3.7 


3.8 


4.1 


5.5 


12.000 


1.6 


8.3 


4.2 


4.4 


4.8 


6.5 


16.000 


1.7 


6.8 


5.3 


5.7 


6.2 


8.5 



TABLE XVII. Here we show the constituent lattice and continuum pieces of 1?^^ in the regions noted in bold in table XVI 
It can be seen that the cancellation between the lattice and continuum one-loop correction is resulting in unstable behaviour 
for the overall scale. 



rrib 


rUc 


aqlAT 


^Qcont 


Zlat 


^cont 


aqlii 


2.0 


1.0 


1.73 


0.74 


0.3039 


-0.1869 


6.8 


2.0 


2.0 


1.63 


1.04 


0.1453 


-0.2122 


0.4 


3.0 


1.2 


1.66 


0.99 


0.1992 


-0.1683 


27.9 


3.0 


1.5 


1.62 


1.10 


0.1501 


-0.1869 


0.2 


4.0 


1.2 


1.63 


1.14 


0.1795 


-0.1371 


5.1 


4.0 


1.6 


1.58 


1.32 


0.1200 


-0.1683 


0.8 


6.0 


1.8 


1.50 


1.71 


0.0878 


-0.1371 


2.1 


8.0 


1.6 


1.48 


1.86 


0.0971 


-0.0804 


0.5 


10.0 


2.0 


1.40 


2.33 


0.0660 


-0.0804 


24.0 
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